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Abstract 
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| We evaluate the radion and Higgs masses in the gauge-Higgs unification models on 

the warped geometry, in which the modulus is stabilized by the Casimir energy. We 
analyze the one-loop effective potential and clarify the dependences of those masses 
on the Wilson line phase #h- The radion mass varies 1-30 GeV for 0.06 < sin#H < 0.3, 
while the Higgs mass is 150-200 GeV and depends on On only logarithmically. The 
radion couplings to the standard model particles are sensitive to the warp factor, and 
$_i ' are too small to detect at colliders in the region where the five-dimensional description 



is valid. 
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1 Introduction 



The gauge-Higgs unification scenario is an interesting candidate for the physics beyond the 
standard model, which was originally proposed in Refs. [TJ |2] and revived by Refs. [3], H] as 
a solution to the naturalness problem. In this class of models, the Higgs mass is protected 
against large radiative corrections thanks to a higher- dimensional gauge symmetry [5]. 
The models are characterized by the Wilson line phase along the extra dimension 8 ft. The 
electroweak symmetry breaking occurs when sin^H 7^ or sin ^ 0, depending on the 
models. The fluctuation around the vacuum expectation value (VEV) of 9ft corresponds 
to the physical Higgs boson in the standard model. 

This scenario has been first investigated in the flat spacetime [6j [7], and extended to 
the Randall-Sundrum warped spacetime j8]. The models in the latter can solve some 
problems that exist in the former case. The masses of the Higgs and the Kaluza-Klein 
(KK) modes are enhanced by a logarithm of the large warp factor [9J so that they can 
evade the experimental lower bounds, and the large top quark mass can easily be realized 
only by the localization of the mode functions in the extra dimension [TU]. Furthermore, 
such models have phenomenologically interesting features [10]- [IS]- Hence, we will focus 
on the Randall-Sundrum spacetime as a background geometry in this paper. 

When we work in extra-dimensional models, the stabilization mechanism for the size of 
the extra dimension, which is often called the modulus or the radion, must be considered. 
One of the simplest mechanisms for the modulus stabilization is proposed in Ref. [T7]. A 
five-dimensional (5D) bulk scalar field plays an essential role for the stabilization in this 
mechanism. The modulus can also be stabilized by the Casimir energy of the bulk fields. 
This possibility has been discussed in many papers [TS]- [21], and it has been shown that 
5D gauge and fermion fields that spread over the bulk are essential for the modulus stabi- 
lization [22] ■ Thus the latter mechanism is more economical in the gauge-Higgs unification 
scenario because the bulk gauge and fermion fields already exist in the theory and no extra 
bulk scalar fields need not be introduced just for the stabilization. 

In our previous work [23], we discussed the modulus stabilization by the Casimir energy 
in the model proposed in Ref. [2] , in which the Wilson line phase is dynamically determined 
as #h = f • We found there that the brane kinetic terms for the gauge fields are necessary 
for the modulus stabilization, and the radion mass is 0(1 GeV). Although this model 
has phenomenologically interesting features [IS], the electroweak precision measurements 
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disfavor #h = f according to the analysis in Ref. [10J. Besides, it is a nontrivial task to 
clarify the #H-dependence of the radion and Higgs masses because the effective potential V e s 
depends on parameters that control the VEV of #h i n a complicated way. Therefore, in 
this paper, we will extend our previous work [23] to the case that #h can take small values 
and clarify the 6*h- dependence of the radion and Higgs masses by evaluating V e s. We will 
also discuss the experimental constraints on the radion mass. 

The paper is organized as follows. In the next section, we provide a brief review of the 
model in Ref. [H] focusing on the matter sector, and show the one-loop effective potential 
for the radion and Higgs fields. In Sec. [3j we extend the matter sector of the model to 
realize small values of 9n, and see how the effective potential is modified by such extensions. 
In Sec. HI we estimate the radion and Higgs masses as functions of #h, and comment on the 
experimental constraints on the radion mass. Sec. |5]is devoted to the summary. We define 
some functions useful for our analysis in Appendix |A] show an approximate form of the 
effective potential in Appendix (Bj and provide some useful expressions for the numerical 
calculation in Appendix 

2 50(5) X U(l) x model 

In this paper, we consider the gauge-Higgs unification models based on a 5D 50(5) x U(l)x 
gauge theory. This class of models was first discussed in Ref. [10], and several similar models 
with different matter sectors have been studied so far (TTJ [T5J |2S] . In our previous work |23j . 
we considered a model proposed in Ref. [H] as the simplest example. We start with a brief 
review of this model, focusing on the matter sector, and extend it later. 

We assume the 5D warped spacetime compactified on an orbifold S 1 jZi [8] as a back- 
ground geometry. The background metric is given by 

ds 2 = G MN dx M dx N = e~ 2a ^r]^dx^dx v + dy 2 , (2.1) 

where M, N = 0,1,2,3,4 are 5D indices and 77^ = diag(— 1, 1, 1, 1). The fundamental 
region of S 1 /Z 2 is < y < L. The function e a ^ is a warp factor, and a(y) = ky in 
the fundamental region, where k is the inverse AdS curvature radius. The orbifold has 
two fixed points y = and y = L, which are called the UV and IR branes, respectively. 
The gauge symmetry is broken to SU(2) L x SU(2) R x U(l)x at the IR brane, and to 
SU(2)l x U(1)y at the UV brane by boundary conditions [TU]. In order to stabilize the 
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modulus, we need brane-localized kinetic terms for the gauge fields [22]. Thus we introduce 
the following terms on the IR brane. 



/•kin _ o /~r 



4k 



tr {F (C) F (G^ } _ ^ tr {F (A) Fi A^ } _ £ F W F Wr] S (y - L) 



4k' 



(2.2) 

where \f—g = det(g^), g^ v is the 4D induced metric on the IR brane, F^^, F^ Nl F^f^ are 
field strengths for the SU(3)c, 5*0(5), U(l)x gauge fields, and k c , k w , k x are dimensionless 
constants. For simplicity, we do not consider kinetic terms on the UV brane nor brane 
kinetic terms for the 5D fermions, and assume that the following. 



2.1 Matter sector 

We introduce 5D fermions (i = 1,2, •■■) belonging to the vectorial representation of 
50(5) as matter fields. The 5D Lagrangian in this sector is given by 



-G 



{i^r^** - iM m e{y)*i%} + 



(2.3) 



where G = det(GMTv), F N are 5D gamma matrices contracted by the funfbein, T>n is the 
covariant derivative, and e(y) is a periodic step function. The ellipsis denotes the gauge 
sector. 

It is useful to express the 50(5) vector \I/ = (ipi, ■ ■ • , 5 )* as 



where 







■011 i> 
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,021 
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(2.4) 



(2.5) 



y 4>21 4>22 J V% 

is a bidoublet and -05 is a singlet for 5?7(2) L x 5£7(2) R . For example, the third generation 
of the quark sector comes from two multiplets and \&2, which are expressed as 







Qi = 




Q2 = 







q 



t' 



x 



b' 



(2.6) 



The orbifold parities for them are listed in Table HI The subscript R denotes the 4D 
right-handed chirality defined by 75 = T 4 . The left-handed components have the opposite 
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QlR 


QR 


t'n 


Q2R 


Q3R 


V R 


(-,-) 


(- -) 


(+,+) 


{--) 


(- -) 


(+,+) 



Table I: The orbifold parities for the fermions. The left and the right signs in the parenthesis 
denote parities around y = and y = L, respectively. The left-handed components have 
the opposite parities to the right-handed ones. 

parities to the right-handed ones. On the UV brane, we can introduce brane-localized chiral 
fermion fields and change the boundary conditions there, just like we did in Ref. [14J. The 
resulting boundary conditions on the UV brane are 

Qil = c^qi + s^Q 2 l = Q3L = 0, 

-s w q R + c w Q 2 r = 0, t' L = b' L = 0, (2.7) 

at y = 0. Here = smu, c u = cosu, and u is a mixing angle determined by the ratio of 
the boundary mass parameters. (Since qi and Q%l have the same quantum numbers for 
<S77(2)l x U(l)y, they can mix on the UV brane.) 



2.2 Mass spectrum 

The mass spectrum {m n } in the 4D effective theory is determined as solutions to the 
equation, 

P / (A„,^ H )=0, (2.8) 

where I specifies the sectors, and X n = m n /k. The functions //(A; #h) are written in terms 
of the Bessel functions, and listed in Appendix A of Ref. [23] E (See ( 1B.1I) . for example.) 

The W and Z boson masses are obtained as the lightest solutions to p w (Xw', Or) = 
and p z (Xz] On) = 0, and are approximately expressed as 

ke~ kL sin 0u 



m w = kX w 



y/kL + K 



1 + s 2 

m z = kX z ~\l ^ke~ kL sin R . (2.9) 

kL + K 

Here, = Qb/ \fg\ + 9%, where g& and gs are the 5D gauge couplings for the SO (5) and 
U(l)x gauge fields. 



1 They are denoted as pi(X) in the notation of Ref. [53], and we need to generalize the expressions there 
to the case that M*i ^ M^i for the quark sector. 
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2.3 Effective potential 

The one-loop effective potential V e s for the radion and Higgs fields is calculated by the 
technique in Ref. |19] . As mentioned in Ref. [22], only the fields that spread over the bulk 
can give sizable contributions to V e R. In our model, such fields are the gauge fields and the 
quark multiplets in the third generation^ 

Now we promote dimensionless constants kL and #h to 4D dynamical fields p(x) and 
6n(x). Then V e g is expressed as the following form. 

V eS (cp, H ) = kv + e-^r m + e~^V{^ H )1 , (2-10) 

where 

Here rjj = (|) for bosons (fermions), iVj is a number of degrees of freedom for a 
particle in sector /. The functions fC 1 (w) and Z 7 (u>) are expressed by products of the 
modified Bessel functions e~ tcm K a (w) and e ll3n Ig(w) respectively^ and defined so that 
p l {iwe-*] 0)/IC I {we-^)l I {w) becomes a product of {l - e iM)7r ^^3p^}- ( See Ap- 
pendix B of Ref. [23].) The dimensionless constants ttjv and tir, cannot be determined in 
the context of the 5D field theory^ 

In our model, V e g is approximately expressed as 

V eS {<p, 9 B ) ~ Ko(^) + V 2 (<p) cos 2 H , (2.12) 

where Vo((p) and V^y?) are independent of #h- From the stationary condition for 9n, we 
obtain 

sin2# H = 0. (2.13) 

Thus we have #h = 0, ±|, it as candidates for the vacuum. As mentioned in Ref. [H], the 
contribution from the gauge sector prefer the symmetric phase #h — 0,7r, while that from 
the fermion sector does the broken phase 6*h = ±f • In fact, due to a large contribution 
from the top quark sector, #h = ±f is selected as a vacuum, and the electroweak symmetry 
is broken. 
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A fermion field spreads over the bulk when its bulk mass is close to k/2. 



3 The definition of Kp{w) and Ip(w) are given in ()A.4|) . 

4 The divergent one-loop contributions to them are given in (B.fO) of Ref. [23]. They are absorbed in 
the renormalization of the tensions of the UV and IR branes, respectively. 
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Qr 


t R 


b R 


(+,-) 


(-,+) 


(-,+) 



Table II: The orbifold parities for the components of ^3. The left-handed components have 
the opposite parities to the right-handed ones. 

3 Extensions of the model 

According to the analysis of Ref. [10], the VEV of On must be small, i.e., 

sin H < 0.3-0.5, (3.1) 

from the constraint on the oblique parameter S < 0.3. 

To realize a small value of On, the matter sector has to be extended. The simplest 
extension is to introduce an additional fermion multiplet ^3 that belongs to the spinorial 
representation of SO (5) and whose U(l)x charge is 1/6. It is decomposed as 



#5 



Q, t, b 



(3.2) 



where Q, i and b transform as 2i/ 6 , l 2 /3 and I-1/3 under S'Z7(2)l X U(l)y- The orbifold 
parity of each component is assumed as shown in Table [III This sector consists of the 



Q 



EM 



and Qem = — k sectors, and their mass spectra are determined by (12. 8p with 



P 



^(XM = „*<-»">(A;fa) = ^j^W^-jW " (3-3) 



where the function F£J\) is defined by (lA.lft . and c 3 = M^/k. This can easily be 
obtained in the usual procedure to determine the mass spectra in the warped spacetime 
(see, for instance, Ref. [12]). Note that the period of the spectrum is 2tt, which is twice of 
those in the other sectors. If the bulk mass is close to k/2, a contribution from \&3 to 
V e fi is sizablejf] The approximate expression of V e s is modified as 



V eS {v, Oh) - V (v) + V 1 {(p) cos^h + V 2 {ip) cos 2 9 H . 



Now we have a linear term for cos^hIj Then, we find a new stationary point of V e s, 



cos 9u 



Vi 



W2 



(3.4) 



(3.5) 



5 For the anomaly cancellation, additional fermion multiplets are required. However, they are irrelevant 
to the current discussion unless their bulk masses are close to k/2. 

6 Vo is also modified from those in (|2.12j) by the ^-contribution. 
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if |Vi| < 2 [ V2 1 - When V^y?) > 0, this gives a global minimum of V e s for a fixed value of if. 
The value of (I3.5P is controlled by the bulk mass Mq 3 . In fact, a small value of #h can be 
realized by choosing it such that M^i ~ M^ 3 . 

There is another extension of the matter sector. In Ref. [11], two additional fermion 
multiplets and ^2 are introduced, which have the same quantum numbers as and 
respectively. The orbifold parity at the UV brane for each component of ^1 and \&2 
are the same as those of and ^2, while the parities at the IR brane for the former 
are opposite to those for the latter. Then the following mass terms are allowed on the IR 
brane. 

££T = 2^g [iCi (QilQir + Mr) + %JrV l 

-NCa (Q2LQ2R + QslQsr) + i&VrVl + h.c. | S(y — L), (3.6) 
where £1,2 and £1^ are dimensionless mass parameters, and 



Qi, q, t' , ^ 2 = Q2, Q3, % . (3.7) 



In general, q^ and Qil can mix with qi and Qil on the UV brane since they have the same 
quantum numbers although it is not considered in Ref. [11] for simplicity. 

The boundary mass terms in (13.61) relate ^1 and ^2 with SE^ and \j/ 2 through the equa- 
tions of motion. This induces a quartic term for cos^h in V e g. Namely, the approximate 
form of V c s now become 

V eS ( V , H ) - V (<p) + V 2 (ip) cos 2 U + Vt(<p) cos 4 U . (3.8) 

In this case, a solution of 

C ° s2 ^~2l| (3 - 9) 
can be a candidate for the vacuum value. In fact, it becomes a global minimum of V e s for 

a fixed value of if when < — V 2 < 2V4 is satisfied. 

We can also extend the matter sector by introducing 5*0(5) tensor multiplets [Til 123] 

with boundary masses among the bulk fields. Also in this case, V e s has the form of (13. 8p . 



4 Modulus stabilization 



4.1 Scalar mass matrix 



From the stationary condition for ip, we obtain 



rm = 9 ~f- - V. (4.1) 
8 



Since tfr cannot be determined within our setup, it should be treated as an input param- 
eter. In order for the 5D description to be valid, the 5D scalar curvature IZ^ = —20k 2 
must satisfy the condition |7£g| < Mf, where M 5 is the 5D Planck mass [26]. Namely, 
£ = M 5 /k > 4.5. For a sufficiently large warp factor, this means that 

k ~ J^M Ph (4.2) 
where Mpi is the 4D Planck mass. By using (I2.9p . we obtain 

e kL VkTT^ ~ — - sin fl H £ 4.5 x 10 15 sin H - (4.3) 



m w V C 3 

This means that the warp factor has an upper bound depending on #h- An extremely 
small value of sin^n does not allow e kL ~ 10 15 . Besides, it requires a fine tuning among 
the model parameters, such as M^j (i = 1,2,3). So we focus on a parameter region where 
0.06 < sin#H < 0.3, and take the allowed maximal value e kL = 5 x 10 13 in the following 
analysis. From ( 14.11) with the explicit expression of V, we can see that this value of the 
warp factor is realized by an 0(1) value of tjr. Here ( 14. II) is evaluated at the value of #h 
determined by ( 13. 5 p or (13. 9p . The other constant Tuv is determined by the condition that 
the cosmological constant in the 4D effective theory vanishes. 

The bulk masses for the fermions control the profiles of the zero-modes and their mass 
spectrum in the 4D effective theory. In the absense of the boundary mixing, the zero-mode 
mass eigenvalue is larger (smaller) than when c < ~ (c > ~), where c is a ratio of the 
bulk mass to k. (See Fig. 2 and Table 1 in Ref. [12] •) The situation is more complicated 
in our case because of the boundary mixing parametrized by in (12. 7p . For example, in 
Ref. [H], M^i = Mq,2 = 0.43/c and = 1 — |(mb/m t ) 2 are chosen to reproduce the top 
and bottom quark masses. In this paper, we will choose M^i = 0.43/c, M*2 = 0.53/c and 
3^ = 0.86 as an example. In the case of ( 13 .4p . M^, 3 determines the value of 6*h- 

Once the warp factor is given, we can discuss the stability of the vacuum. By using the 
stationary conditions, the second derivatives of V e g at the minimum are given by 

^4 e -4fcL 



1 





L4 -4fcL 





(4.4) 
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where the symbol |o indicates that the quantity is evaluated at the minimum of V e g. These 
provide the mass matrix for the fluctuations (p = (p — kL and 6-^ = 6u — (On). Note that 
we have to canonically normalize these fluctuations in order to discuss the physical masses. 
The canonical normalization for them are given by 



, 3M| . 



4k 



Then the (squared) mass matrix for r and h is calculated as 

2 2 
2 / m rr m rh 



Scalar = 7 I I 



where 



2 2 
m rh m hh 



h(p 2kL — 1 \ l.5 p -2kL . 

m 2 = ^ -!-d 2 Va\ - — (d 2 V-4dV 

n 2kL 1 \ „ L,A -2kL 



m 



g A (e 2kL - 1) a a T , . , 



rh 2v^M| * '° 327r 2 v/3M| * 



o 







2/2fcL_i\ 2l3 -2fcL 

2 _ ffM e 1 J o2 t/ I ~ ^ e t A 



(4.7) 



4.2 ^-dependence of various mass scales 



Now we express each parameter in terms of the 4D ones, i.e., Mpi, my/ and the 4D S77(2)l 
gauge coupling g 4 . First we should note that k is a function of #h through (12.91) since we 
take m w as an input parameter. Thus M 5 and g A also depend on # H as 



M 5 = i 2fc 2fc L M Pi - 2e VfcL + K Mg 1 m w 
1 — e sm 

1/2 



gWkL + k ( \fkL~+K sin H \ 
c/a = 7= ^ 5-4 ^ • 4.8 

Let us first consider the case of (13.41) . As shown in Appendix[Bl V"(^, #h) has the following 
approximate form for tp ^> 1. 

V^, H ) = W Un + M cos™ H , (4.9) 
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where u n and v n are constants, and uq = 0(— 5), u% = u% — 0(3), t>o = t>2 = 0(10) and 
t>i = 0. Then, by using (14. 8 p and (14. 9ft , the expressions in ( 14. 7ft become 

„2fcL„ 



2 ^ ' "'ir , 2 



(v + V 2 COS 2 H ) 



2 ^ (J\< U n> : \\- 



m rh ^ — ;= — — hr^2Cos0 H) 

^ 2 ^#^4^ 2 + ^y (410) 



24ir 2 M£ } sin 4 # H 

16V6tt 2 M P i sin 2 # H 
glkL 2 m 2 w / v 

* h - 32vr2 r 2 + fcz; • 

Here we have assumed that k <C fc-L ~ 32. Since e fcL miy/Afpi = 1.7 x 10~ 3 <C 1 for 
our choice of the warp factor, the mixing angle between the radion and the Higgs boson 
is negligible. Notice that there is no kinetic mixing between them in our model, which 
originates from the curvature-scalar mixing term y/— g^TZ^ (g)H^H5(y — L) [27] 128] (£ is a 
dimensionless parameter, TZ^g) is the 4D Ricci scalar for the induced metric g^ v and H is 
the Higgs doublet), because such a term is prohibited by the 5D gauge symmetry when H 
is a part of the 5D gauge field. 

Thus the radion and Higgs masses m rac i and run are roughly estimated as 



m w 

m rad ~ z—7=—r: — . 9 ~ Vvo + v 2 



^/0(2O) x 1(T 2 GeV 



2 v / 6vrM P1 sin 2 ^H sin 2 6 R 



^ + ^~^)x Ma , (4, i; 



We have used that cos6>h ~ 1. 

We obtain a similar result also in the case of ( 13.81) . Now V(f, #h) is approximated as 

V(<p, H ) = V ( u 2n + ^) cos 2 " 9 H , (4.12) 

where uq = 0(— 5), u 2 = 0(3), = 0(1), t> o = v 2 = 0(10), and v& = 0. Then, m ra( j and 
mi are estimated as 



,fcL 



n? 



v/0(2O) x 10" 2 GeV 



2v / 6vrM P1 sin 2 ^ H sin 2 9 R 

mH „ — — — -\/u4 cos ~ x 180 GeV. (4.13) 

2v 2vr 

The radion-Higgs mixing is negligible also in this case. 
The typical KK mass scale ttikk is estimated as 

k VkLm w 450 GeV , , , 

mKK = ~ki — T —q — - —q — • 414 
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Figure 1: The masses of the radion and the Higgs boson as functions of #h- The solid, 
dotted, dashed lines represent the case of k = 5.0, 3.0, 1.0, respectively. 
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Figure 2: The mass of the first KK mode for ^3 as a function of #h- 

Now we show some numerical results. As a specific example, we consider the case of 
(13 .4p . Fig. [TJ shows the radion and Higgs masses as functions of #h- The lightest KK mode 
comes from \l/3 when k < 20. Fig. [2] shows its mass as a function of #h. From these plots, 
we can read off the ^-dependence of each mass as 

0.04 GeV 



m Tad 



sin 2 #h 



m H ^ (150 - 10 In sin H ) GeV, 



(i) 



k 



230 GeV 



(4.15) 



2e kL sin #h ' 

for k = 1.0. The radion-Higgs mixing angle is less than O(10~ 4 ) for sin^n > 0.06. These 
results are consistent with the rough estimations (14. lip and ( 14. 14j) . 

Here we comment on the strength of the brane kinetic terms. We find that the modulus 
stabilization requires k > 0.5. For such values of k, the Higgs mass is larger than 175 GeV 
at #h = 0.06. Therefore, within the parameter space we consider, there is a region that is 
consistent with the latest excluded region 158 GeV< m H <175 GeV by the Tevatron [30] . 
Since the zero-modes for the gauge fields have (at least approximately) flat profiles, the 
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4D and 5D gauge couplings g± and g$ are related through 



2 v ' ' ~TT ) ' ( 4 -!6) 



ft 

Thus the dominant contributions to g± come from the bulk terms and the brane kinetic 
terms only give small corrections for an 0(1) value of ft. On the other hand, the brane 
kinetic terms with 0(1) ft affect the masses of the KK modes and their couplings to the 
zero- modes in a sizable way, and their impacts are discussed in Ref. [23j [29] , for example. 

4.3 Experimental constraints on the radion mass 

Finally we comment on the experimental constraints on the radion mass. The couplings 
of the radion to other fields are obtained by expanding the 4D effective action in terms 
of <p = r/A r , where A r = ^/3Ml /(k(e 2kL - 1)) ~ y/Ee~ kL Mpi. In the original Randall- 
Sundrum setup, in which all the standard model particles reside on the IR brane, they are 
expressed as [3TJ [32] 

C eS = -^-rT^ + ..., (4.17) 

where T^ v is the energy-momentum tensor of the standard model. Namely, the radion 
couples to particles on the IR brane just like the standard model Higgs does with an extra 
factor v/A r , where v = 246 GeV. To particles propagating in the bulk, the radion couplings 
deviate from (I4.17p . but are of the same order of magnitude [33J ES]- Hence the radion 
couplings are very weak because v/A r — 2.1 x 10~ 3 for our choice of the warp factor. Thus 
the radion mass is not constrained from the collider experiments^ It can also be bounded 
from below by the consideration of the neutrino oscillation inside the supernova [SB], but 
this lower bound is much lower than 1 GeV in our case. 

In addition to the above tree-level couplings, the radion couplings to the photons and 
to the gluons also receive sizable one-loop contributions^ These couplings are important 
for the production and the decay processes of the radion. They can be enhanced compared 
to the corresponding Higgs couplings [^| times v/A r . (See, for example, Ref. [H].) However 



7 The radion with a mass in the range: 12 GeV < m ra d < 90 GeV is excluded if (v/A r ) 2 > 0.01 [36 , 
and it is excluded for m rad < 12 GeV if (v/K r ) 2 > 0.1 [37] , 

8 In contrast to the original Randall-Sundrum model, there are tree-level contributions to these couplings 
in our models because the gauge fields propagate in the bulk. 

9 The Higgs couplings to the massless gauge bosons are discussed in the context of the gauge-Higgs 
unification in Refs. l39l 1401. 



13 



such enhancements are insufficient for compensating the suppression factor v/A r , and for 
discovering the radion at the Large Hadron Collider. 

The situation does not change so much even if a larger value of e kL is chosen. From the 
requirement that M 5 /k > 4.5 and sin^n < 0.3, the warp factor is bounded through ( 14. 3 p 
as e kL < 2.3 x 10 14 for an 0(1) value of k. Thus the suppression factor v/A r cannot be 
larger than 0.01, which is still too small to detect the radion at the colliders. 

In the case that the fermion mass hierarchy is realized by the wavefunction localiza- 
tion in the extra dimension, the experimental bounds on the flavor- changing processes can 
provide stronger constraints on m ra d and A r . According to the analysis of Ref. [35J, they 
should satisfy A r m ra d > 2.3a^ s TeV 2 , where dd s is a dimensionless constant that parame- 
terizes the flavor violatio and its typical values range between 0.03 and 0.12. Thus the 
most stringent bound on m ra d is that m ra d > 2.7 GeV for e kL = 5 x 10 13 and m ra d > 10 GeV 
for e kL = 2.3 x 10 14 . 



5 Summary 

We have considered the modulus stabilization in the gauge-Higgs unification scenario, and 
estimate the radion and Higgs masses. Through the ^-dependences of M 5 and g& in 
( 14. 8p . various mass eigenvalues depend on #h nontrivially. We found that the masses of 
the radion, the Higgs boson and the KK modes all have different ^-dependences. In order 
to see them explicitly, we considered two classes of models which correspond to different 
extensions of the model in Ref. [2] . Qualitatively, we have the same results in both classes. 

10- 2 GeV 



^rad 



y/0{2G) 



sin 2 6*h 
m H ~ ^0(3) x 90 GeV, 
0(500 GeV) 

»KK ~ — n , (5.1) 

sin R 



and the radion-Higgs mixing is negligible. As mentioned in our previous work [23], the 
boundary kinetic terms for the gauge fields are necessary for the modulus stabilization. 
An 0(1) value of k in ( 12. 2ft is enough to stabilize the modulus. In contrast to the model in 
Ref. [2] , the Higgs couplings to other particles do not deviate very much from the standard 
model values when #h <C 1 • 
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It is determined by the bulk mass parameters for the quark multiplets in the first two generations. 
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The radion couplings to the standard model particles are suppressed by factors of 
0(v/A r ) compared to the corresponding Higgs couplings. Within the allowed region of the 
parameter space, such suppression factors are less than 0.01. Therefore collider experiments 
do not impose any constraints on m ra d. However the experimental bounds on the flavor- 
changing processes can provide stronger bounds on m ra d if the fermion mass hierarchy 
stems from the wavefunction localization. Such bounds narrow the allowed range of #h in 
some cases. 

Cosmological impacts of the radion physics is an intriguing subject. For this direction, 
we might need to extend the works by Refs. [12| H3j HI] to deal with the one-loop effective 
potential at finite temperature in the Randall-Sundrum background. This is one of our 
future projects. 
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A Definitions of functions 

Here we define functions that are useful to express //(A; 6^) in ( 12. 8 1) and the effective 
potential. First we define the following functions from the Bessel functions. 

F^(X) = J a {\)Y${\z L ) - Y a (\)J$(\z L ), (A.l) 

where z L = e kL , and 

J%{u) = J p {u) - nuJ p+l {u), Y^{u) = Y p (u) - kuY p+1 (u). (A.2) 
For calculations of the effective potential, we also define 

KA W ) = I a {w)K*{wz L ) - e~< a -®*K a (w)I%(wz L ), (A.3) 

where 

I$(u) = I p (u) + KuIf, +1 (u), K K p {u) = Kp{u)-KuK p+1 {u). (A.4) 
Then the following relation holds. 

F^{iw) = --e^-^ 2 F^(w). (A.5) 

7T 

The asymptotic behavior of F£ p(w) for Hew ^> 1 is 

w(z L -l) 

F^(w) = -__ e <0-«>* (1 + nwz L ) {l + O (w' 1 ) } . (A.6) 
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B Approximate form of effective potential 



Here we show the approximate forms of (14. 9p . As mentioned in Ref. [22], the dominant 
contributions to the effective potential come from the gauge fields and the fermion fields 
with bulk masses that are close to k/2. 

As an example, let us consider a sector whose mass spectrum is determined by 

/(A; 9 U ) = F^ W (A)F Q ° Q (A) - = 0, (B.l) 

where a is close to one. In fact, a = 1 for the gauge sector, and a = M^/k + 1 for the 
fermion sector. Then, the integrand of (12.111) can be approximated for ip ^> 1 as 

ln ptjiwe-^-en) = p J (zwe-^0) p 1 \iwe^- H ) 
n /C / (we-^)X / (w) n /C / (we^)X / H n ^(iwe-v-Q) 



lf a _i(«;e-' 6 )/*_ 1 («;) J I ]f a (ice-f)J tt (ifi) 

2w 2 F a K _ lia _ 1 (toe^)F a a («;e^) 
7 a _! (iwe-v) Kl_ x (w) sin 2 H 



K a -.i(we-v) /;_jH 2w 2 F Q K _ liQ _ 1 (u;e-^)F Q Q (^e-^) 



r(a-i)r(a) \ /S-iM ^-i(«V«H J ' ^ j 
I (-W, sin 2 g H \ ( a = l) (B 2) 

sin^ H | (a<1) 

where kL in the definition of //(A; #h) is replaced by tp, and r(a) is the Gamma function. 
We have used that \a — 1| <C 1, and assumed that w = 0(1) because the above func- 
tion exponentially decays for w 1. Thus the contribution of this sector to V((p, Or) is 
negligible when a > 1. When a < 1, it is estimated as 

Jo /C J (we tp )l 1 (w) 



^ELo^cos^h, (a = l) 

(-) 2% ^ELo^cos^h, (a<l) 



(B.3) 



where 



v = N! dw W 3 < — + 







J^M 2wI^ 1 (w)I a (w) 



w 2 



0, ^ = —Ni / — . - . . . (B.4) 

1 2/-_ 1 (w)/ Q H 1 ; 



n 
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Here we have used that <p ^> |ln(w/2) + 7I in the integration region that gives dominant 
contributions. For k — 1.0 and a = 1, for example, Vq ~ 8.6 and v\ ~ —5.4. In general, we 
can see that v\ = O(Nj) and v\ = 0(—Nj)\^j For the gluon and photon sectors, in which 
p^ =< ^'^(A; $h) = A.F K (A), tlie terms proportional to sin 2 in flB.2j) are absent. Thus we 
find that v* = -0{Ni) and v[ = v{ = 0, where N G = 24 and iV 7 = 3. 

For the ^-sector, in which p 7 (A; 9n) is given by (13 . 3 p . we can estimate V I=iS,3 (tp, #h) in 
a similar way and find that Vq s = <D(N* 3 ) = (9(24), v* 3 = (9(-24) and v% = 0. 

As a result, we obtain 

1 

2sin(7rai) . w Wl 2 s 2sin(7ra 3 ) , * 3 ^ 3 x 

~ ( W + v 2 COS #H ) ( Vq + Vi COS 6» H J 

7T 7T 

7=W,Z ^ i=G,7 ^ 

= E f u ™ + ~) cos "^h, (B.5) 
where uq = 0{— 5), «i = M2 = (9(3), t> = v 2 = (9(10) and V\ = 0. 

C Expressions for numerical calculations 

Although the approximate expressions of the mass eigenvalues in (14. lip or (14.131) is useful 
for the order estimation of the mass eigenvalues, we need more accurate expression for the 
numerical calculation. 

Here we consider the case of (13.41) as an example. Then, since the functions //(A;#h) 
have the form of 

p 7 (A; H ) = Po(A) + p{(A) cos^h + pi{X) cos 2 H , (CI) 
we can expand the integrand of (12.111) around cos #h = c (c is a constant) as 

In// = In {po + p[c + P2C0} - cj[cq + a\<? Q 

+ [a[ - 2<4cq) cos 6 n + o[ cos 2 6 U + O ((cos 6 U - c ) 3 ) , (C.2) 



11 For the fermion sector, the forms of p I (X;9-a) are more complicated than (|B.1[) . but the above rough 
estimate does not change much. 
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where 

^ = p{ + 2p I 2 c ^ _ 2p I pj - {p[f - 2p{p I 2 c - 2 (pj) 2 c\ 

Gl ~ pi + p{c + p[cl ' ° 2 - 2{pi + p[c + pjc 2 Q f 

Thus V(ip,9n) in (I2.1ip is expanded as 

V((p, H ) = V Q ((p; c o) + V x {(p\ c ) cos 6 U + %{?\ c ) cos 2 6 n + O ((cos# H - c ) 3 ) , (C.4) 

where 

V (<p; c = >,(-) m Ni / dww 6 < In ^ -j- — 

Jo [ 1C l (we <^)Z 7 (w) 

- (<t[cq - cr^Cg) {iwe" 9 )} , 

POO 

Mr, Co) = y^H^'Nj / w 3 (a[ - 2a{c ) (iwe^), 
j Jo 

POO 

V 2 (v;c ) = yV-^'iVi / dtoA'f™^). (C.5) 

If we choose a constant Co as a solution to the equation, 

Vi(kL;c ) 

c ° = ~{F7Tt — V (C ' 6) 

V 2 {kL;c ) 

the correction term (^((cos^h — Co) 3 ) in (IC.4I) can be neglected in the calculations of m ra d 
and mH. 
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